Electronic structures of (In,Ga) As/GaAs quantum dot molecules made of dots with 

dissimilar sizes 



O 

o 

X> 

PL, 

l> 
(N 



Lixin He 1 and Alex Zunger 2 

1 Key Laboratory of Quantum Information, University of Science and Technology of China, 
Hefei, Anhui 230026, Peoples Republic of China 
^National Renewable Energy Laboratory, Golden, Colorado 804-01, USA 
(Dated: February 6, 2008) 

Using single-particle pseudopotential and many-particle configuration interaction methods, we 
compare various physical quantities of (In,Ga)As/GaAs quantum dot molecules (QDMs) made of 
dissimilar dots (heteropolar QDMs) with QDMs made of identical dots (homopolar QDMs). The 
calculations show that the electronic structures of hetero-QDMs and homo-QDMs differ signifi- 
cantly at large inter-dot distance. In particular (i) Unlike those of homo-QDMs, the single-particle 
molecular orbitals of hetero-QDMs convert to dot localized orbitals at large inter-dot distance, (ii) 
Consequently, in a hetero-QMD the bonding-antibonding splitting of molecular orbitals at large 
inter-dot distance is significantly larger than the electron hopping energy whereas for homo-QDM, 
the bonding-antibonding splitting is very similar to the hopping energy, (iii) The asymmetry of the 
QDM increases significantly the double occupation for the two-electron ground states, and therefore 
lowers the degree of entanglement of the two electrons. 

PACS numbers: 73.22.GK, 03.67.Mn, 85.35.-p 
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I. INTRODUCTION 

Vertically coupled quantum dotsi^ obtained via epi- 
taxial growth provide a potential scheme for scalable 
nano-structures for quantum computing. In this scheme, 
two coupled quantum dots are used as a basic logic gate, 
via the entanglement of one exciton^ or two electronic 
spins. 4 This proposal for gate operations, requires knowl- 
edge of the detailed physical properties of the "quantum 
gate" made of two quantum dots. Significant progress 
has been recently made^i using quantum dot molecules 
made of very large (~ 500 - 1000 A) electrostatically con- 
fined dots. The limit of large quantum confinement, how- 
ever, requires working with (200 x 30 A) self-assembled 
QDMs. So far, most experiments on self-assembled QDM 
are optical^ and most theories are based on continuum 
models, such as effective mass approximations. 3 These 
simple models ignore or drastically simplify important 
real material properties such as strain, atomistic sym- 
metries and crystal structural effects, band coupling etc. 
Recent studies 7 - show that simplification of such impor- 
tant effects may lead to qualitative changes in fundamen- 
tal physics of the QDMs. 

Previously, we have studied homopolar QDMs made of 
two identical quantum dots, using single-particle pseu- 
dopotential method and many-particle configuration in- 
teraction methodi££ We have studied electron localiza- 
tion, double occupation rate and two-electron entangle- 
ment using a new formula for measuring the degree of 
entanglement formula for two indistinguishable fermions. 
We found that even geometrically identical dots in the 
QDMs lead to electronic asymmetry due to the strain ef- 
fects. However, experimentally it is hard to control the 
shape, size and compositions of individual dots within 
the QDMs, so in practice, the QDMs are never made of 
identical dots. Actually, the top dots are tend to be larger 



than the bottom dots due to the strain effects*^ Indeed, 
the measured difference in exciton energy due to dot-size 
difference is about 4 meV— for two vertically coupled 
dots that are 20 nm apart. Sometimes, the two dots 
are intentionally grown different so that they can be ad- 
dressed separately^ To provide quantitative comparison 
to experiments, considering the effects that asymmetry of 
quantum dots within the molecule, we studied the QDMs 
made of (In, Ga) As/GaAs quantum dots of different sizes 
(hetero-polar QDM). 

In this paper, we study systematically the elec- 
tronic properties of hetero-QDMs, including their single- 
particle molecular orbitals, many-particle states, double 
occupation and entanglement of two-electrons, and com- 
pare them to those of homo-QDMs. We found that while 
at short inter-dot distance, the electronic properties of 
hetero-QDM and homo-QDM are similar, they differ sig- 
nificantly at large inter-dot distance. This difference may 
have substantial impact in implementation of quantum 
gates. 



II. METHODS 

Figure Q] shows the geometry of a hetero-QDM, con- 
sisting a pair of 3 nm tall InAs dots in the shape of 
truncated cones, grown on two-dimensional InAs wetting 
layers, embedded in a GaAs matrix. The inter-dot sep- 
aration d is defined as the distance between the wetting 
layers of top and bottom dots. We choose the base di- 
ameter of top dots (labeled as a) to be 20 nm, and that 
of the the bottom dots (labeled as (3) to be 19 nm, mim- 
icking to the fact that experimentally the top dots are 
slightly larger than the bottom dotSfi^^ The composi- 
tion of the dots vary from Ino.5Gao.5As at their bases to 
pure InAs at their top, as determined in Ref. H. We de- 
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FIG. 1: The geometry used in this work for quantum dot 
molecules made of dissimilar dots. We denote the (isolated) 
top dot "q" and the (isolated) bottom dot "/?". Each dot 
has the shape of a truncated cone. The inter-dot distance is 
measured from wetting layer to wetting layer. 
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FIG. 2: Left panel: The single-particle energy levels of molec- 
ular orbitals vs. inter-dot distance. Right panel: The electron 
single-particle energy levels of the isolated dots a and /3. 



note the dot molecules made of dissimilar dots a and (3 
as M a p. We also constructed the homo-QDM, consisting 
a pair of quantum dots 7, which have the average sizes, 
and the same alloy compositions of dots a and (3 in the 
heteropolar dot molecule. We denote the homo-QDM as 
M 77 . 

The single-particle energy levels and wavefunctions of 
M a fj and M 77 are obtained by solving the Schrodinger 
equations in a pseudopotential scheme, 



iv 2 + V ps (r) 



ipi{r) = e % tpi(r) , 



(1) 



where the total electron-ion potential V ps (r) is a super- 
position of local, screened atomic pseudopotentials v a (r), 
and a nonlocal spin-orbit potential V so i.e., V ps (r) — 
a v a( r - Rti,o) + V so . The atomic position {R„, a } 
is obtained from minimizing the total bond-bending and 
bond-stretching energy using the Valence Force Field 
(VFF) model. 12 ' 13 The atomistic pseudopotentials v a 
(a=In, Ga, As) are fitted to the physically important 
quantities of bulk InAs and GaAs, including band ener- 
gies, band-offsets, effective masses, deformation poten- 
tials and alloy bowing parameters, etc?^ Because for 
electrons the spin-orbit coupling is extremely small in 
the InAs/GaAs quantum dots, we ignored this effect. In 
general, including the spin-orbit coupling effect will in- 
troduce mixture of different total spin states. Equation 
(JXJ) is solved in the basis of {<fi rn ^ \ of Bloch orbitals 
of band index m and wave vector k of material A (= In As, 
GaAs), strained uniformly to strain *e following Ref. fl5l 
The Hamiltonian of interacting electrons can be writ- 
ten as, 



Ci a is a fermion operator, ipi— °~ui o~ gi 7r M , ~K g are the 
single-particle eigenfunctions of the i-th molecular or- 
bital, and (J, a'—l, 2 are spin indices. The rjy are the 
Coulomb integrals between molecular orbitals ^ , ipj > ^fc 
and ipi, 



1 kl 



droit 



e(r — r')|r — r'| 



(3) 



The ./,,, = and K tj 



Tfj are diagonal Coulomb and 



exchange integrals respectively. The remaining terms are 
called off-diagonal or scattering terms. All Coulomb inte- 
grals arc calculated numerically from atomistic wavefunc- 
tions^ We use a phenomenological, position-dependent 
dielectric function e(r — r') to screen the electron- 
electron interaction.!- The many-particle problem of 
Eq.([2]) is solved via the CI method, by expanding the 
V-electron wavefunction in a set of Slater determinants, 
\® ei ,e 2 ,-,e N ) = ct^U ■ ■ ■ 4jvl $ o) , where c\. creates an 
electron in the state e» . The i/-th many-particle wave- 
function is then the linear combinations of the determi- 
nants, 



A v (e-L,e 2 , ■ ■ ■ ,ejv) |$ ei , 



(4) 



ei,e 2 ,--- ,ejv 



For the two-electron problems, our calculations include 
all possible Slater determinants of six confined molecular 
orbitals. 



III. BASIC ELECTRONIC STRUCTURES AT 
THE SINGLE-PARTICLE LEVEL 



ijkl a, a' 



n^Ul^k^la, (2) 



where ipia-(r) — ^^^^(r) is the field operator, whereas 



A. Double-dot molecular orbitals 

We first show the electronic structure of isolated dots 
a and (3. The single-dot electron s and p levels of dots a 



3 



and (3 are shown on the right panel of Fig. [2] We see that 
the s-p energy spacing of dot a is e(p a ) — e(s a ) =52 meV 
and that of dot j3 is e(pp) — e(sp) =59 meV, compared 
to 54 meV of dot 7 (not shown). The energy level of sp, 
is slightly (~ 6 meV) higher than s a , because dot (3 is 
smaller than dot a and therefore has larger confinement. 
The p levels of all dots have a small energy splitting due 
to the underlying atomistic symmetry, e.g., 5e(p a ) =6 
meV, and Se(pp) =1 meV. We further calculated the fun- 
damental exciton energy of dot a, Ex(a) = 1153 meV, 
and that of dot /3, E x {(3) = 1159 meV. The energy dif- 
ference in exciton of dots a and (3 is about 6 meV, in 
agreement with experiment^. The fundamental exciton 
energy of the "averaged" dot 7 is Ex{"l) = 1156 meV. 

When two dots a and /3 couple, the bonding and anti- 
bonding "molecular orbitals" ensue from the single-dot 
orbitals. The energy levels of molecular orbitals are 
shown on the left panel of Fig. [2j We show the single- 
particle levels of molecular orbitals^ a g , a u originating 
from s orbitals, and tt u , and Tr g originating from p or- 
bitals. The bonding and anti-bonding splitting A CT = 
e(cr u ) — e(<7 9 ) and A T = e(Tr g ) — e(7r u ) increase with the 
decrease of inter-dot distance, because the coupling be- 
tween the top and bottom dots gets stronger. This pic- 
ture is similar to what we obtained for homo-QDMs. 
However, there is an important difference between the 
homo-QDMs M 11 and hetero-QDMs M a p : in the for- 
mer case, the bonding and anti-bonding splitting A a and 
A w decay to almost zero at large inter-dot distance, while 
in the later case, A CT and A v tend to constants (A<j ~ 7 
meV, Ajr ~ 10 meV here), because the molecular orbitals 
gradually convert at large inter-dot distance to single dot 
energy levels, e.g. the a g levels convert to top dot s or- 
bitals, and a u convert to bottom dot s orbitals, therefore 
the energy splitting between the first and second molec- 
ular states at large distances is approximately the energy 
difference between s orbitals of the top and bottom dots, 
i.e., A CT ~ eisp) - e(s a ) ^ for M a0 . 

Figure [2] shows that at inter-dot distance d =10 nm, 
the molecular orbital levels are about 25 meV higher than 
the isolated dot levels, although the direct electronic cou- 
pling between two dots is much smaller than this quan- 
tity. This energy shift results from the long range strain 
effects experienced by one dot due to the presence of the 
second dots. This effect is missed in EMA-type model 
calculations which ignore strain effects. 



B. Single dot- localized orbitals 

The above discussions pertain to the basis of double- 
dot molecular orbitals. An alternative way to study 
QDMs is to use a dot-localized basis. We have 
demonstrated^ that dot-localized orbitals can be a use- 
ful tool to analyze the QDM physics, including the elec- 
tron double occupation, and two-electron entanglement. 

Dot-localized orbitals Xn can be obtained from a uni- 
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FIG. 3: (Color online) Left panel: The energy levels of dot- 
localized orbitals for QDMs M a p (black solid lines) and M 11 
(red solid lines), and denote the s orbitals of the top 
and bottom dots respectively. The molecular orbitals energy 
levels o g and a u (dashed lines) are shown for dot molecules 
Map- Right panel: s levels of isolated dots a, j3 and 7. 



tary rotation of molecular orbitals, i.e., 

N 

Xr,=22^V,i^i' ( 5 ) 
i=l 

where, ipi is the i-th molecular orbital, and U is a unitary 
matrix, i.e., WU = I. We choose the unitary matrices U 
that maximize the total orbital self-Coulomb energy^*^ 
The procedure of finding hi is described in the Appendix 
B of Ref|9j. As we will show below these dot-localized 
orbitals x-n have the advantage of being only weakly de- 
pendant to the inter-dot coupling. This invariance may 
provide simplified pictures for qualitatively understand- 
ing of the QDM physics. 



1. Single-particle energies of dot-localized orbitals 

The single-particle levels of dot-localized orbitals and 
the hopping (or tunneling) term between two dots can be 
obtained from 

e n = (x v \H \Xn) =^2K,i U V;i e * > ( 6 ) 

i 

tvil2 = (Xvi\Ho\Xv2) = zL^Vui U V2,i e i i ( 7 ) 

i 

where, is the single-particle energy of i-th molecu- 
lar orbital and Ho = ^ icr e Q '0] cr '!/w is the single-particle 
Hamiltonian. Figure [3] depicts the single-particle levels 
ex and e# of the dot-localized orbitals of both top and 
bottom dots, for inter-dot distances d in the range from 
4 nm to 10 nm. (Here, we denote the top dot T and the 
bottom dot -B, to distinguish them from isolated dots a, 
(3 and 7). ex and es of M a p are shown in the black 
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FIG. 4: (Color online) The inter-dot hopping energy 2t (solid 
lines) of hetero-QDM M a p and homo-QDM M 77 . We also 
show the bonding-antibonding splitting A CT of M a p and M 77 



solid lines, and those of M 77 are shown in the red solid 
lines, At large d, the energy difference e# — ~ 6 meV 
for M a p, is close to the value of difference e(s@) — e(s a ) 
between s orbitals of isolated dots a, (3. This energy 
difference gets smaller when the two dots move closer, 
because the energy levels of the top dot rise faster than 
those of bottom dots due to the strain asymmetry. For 
the homo-QDMs M 77 , er and es are almost degenerate. 
The small difference (~ 1 meV) between them is due to 
the strain and alloy effects. We also plot in Fig. [3] the 
energies of molecular orbitals a u and a g in dashed lines 
for M a p. As we see, for d >9 nm, the dot-localized state 
eg of M a p is almost identical to the molecular orbital a u , 
while ex merges with a g , indicating at large d, molecular 
orbitals convert to dot-centered orbitals for M a p. 

The quantity 2t measures the coupling strength be- 
tween the top and bottom dots, and directly determines 
the two-electron properties such as singlet-triplet split- 
ting in the QDM. We calculate this hopping energy be- 
tween the s orbitals of top and bottom dots at different 
inter-dot distances for both M a p and M 77 in Fig. [4j (We 
ignore the orbital index "s" to simplify the notation.) 
We find that 2t(M a p) and 2t(M 77 ) are almost identical 
at all inter-dot distance. However, the hopping ener- 
gies calculated here are much larger than we obtained 
for the pure fnAs/GaAs QDM2, because the alloy QDM 
have much smaller energy barrier between two dots than 
pure QDM. fn general, the quantity 2t does not equal 
to the bonding-antibonding splitting A CT = \J P- + At 2 , 
where 5 = ei&r) — e(es), being the energy difference of 
s orbitals of the top and bottom dots. For homo-QDMs, 
where 8/2t -C 1, we have 2t ~ A a as seen in Fig. 2J How- 
ever, for hetero-QDMs, A CT may be significantly differ- 
ent from 2t, especially at large inter-dot distances, where 
S/2t ^> I, also illustrated in Fig. 0J Experimentally, 19 
one usually measures the bonding-antibonding splitting 



FIG. 5: The Coulomb energies of dot-localized orbitals 
of hetero-QDM M a p. Jtt and Jbb are the s orbital 
self-Coulomb energies of top and bottom dots respectively, 
whereas Jtb are the Coulomb energies between s orbitals of 
the top and the bottom dots. 



rather than the hopping 2t. Therefore, to get the hopping 
energy between two dots, one need to know the energy 
difference 8 of two dots. 



2. Coulomb integrals of dot-localized orbitals 

The Coulomb integrals in the dot-localized basis can 
be obtained from Coulomb integrals of molecular orbitals 
as follows, 

= E K ,i K j u ^ ,k u m rg , (8) 

i,j,k,l 

where are the Coulomb integrals in the molecular 
basis. The direct Coulomb integrals Jtt, Jbb and Jtb 
for M a p are shown in Fig. [5l The Coulomb integrals 
Jtt ~ J<m=2I.4 meV and Jbb ~ J/3/3=22.3 meV, are al- 
most constants at all inter-dot distances, suggesting that 
the dot-localized orbitals are approximately unchanged 
for different inter-dot distance d. Jpp > J QQ , as dot (3 is 
smaller than dot a. The inter-dot Coulomb interaction 
Jtb decay slowly as 1/d. The exchange energies (not 
shown) between the top and bottom electrons is orders 
of magnitude smaller than the hopping energy, and there- 
fore can be ignored in practice. For the homo-QDM M 77 , 
we found that on-site Coulomb energies Jtt ~ Jbb , both 
are very close to the average values of Jtt and Jbb of 
M a a. The inter-dot Coulomb energies Jtb of M a p and 
M 77 are also extremely close. 
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Two-electron states from CI calculation 



Singlet-triplet splitting vs. inter-dot distance 
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FIG. 6: (Color online) Two-electron states for (a) hetero- 
QDM M a p and (b) homo-QDM M T7 , including the singlet 
1 Hg, 1 £g 6 ' states and the 3-fold degenerated triplet 

states 3 E„ as well as two 3-fold degenerated triplet states 

1 1 u ■ 



IV. TWO ELECTRONS IN THE DOT 
MOLECULE 

A. Many-body energy states 



The two-electron-in-a-QDM problem is of special inter- 
est, as it is the prototype of quantum gate using QDMs. 4 
We calculate the two-electron energy levels by the con- 
figuration interaction method using Slater determinants 
constructed from confined molecular orbitals a g , a u and 
7r u , 7r g , which give 66 configurations in total. The two- 
electron energies £ and 3 II tl for hetero-QDMs M a p are 
plotted in Fig. [6^a). To compare with homo-QDMs, we 
show the two-electron states of M 77 in Fig. [DJb). The 
energy levels of M a p are similar to those of M 77 , in the 
following way: (i) The order of the CI levels is unchanged, 
particularly the ground states are still the singlet states 
at all inter-dot distance; (ii) The trend of each CI 
level vs. inter-dot distance d is similar to what we ob- 
tained for M 77 . There are also some differences between 
the hetero-QDMs M a p and homo-QDMs M 77 , especially 
at larger inter-dot distances. For example, in the ho- 
mopolar QDMs, the 1 S„ state is almost degenerate with 
1 T,g at large inter-dot distance, while in M a p, 1 Sg is 
about 13 meV higher than 1 T, U at d=l0 nm. At large d, 
1 E U and correspond to the states that two electrons 
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FIG. 7: The singlet-triplet splitting Js-t vs vs. inter-dot 
distance for hetero-QDM M a p (solid line) and homo-QDM 
M~ ll (dashed line). 



localized on the same dots4i2, The energy difference be- 
tween and 1 E U is due to the size difference of dots 
a and [3. 

The singlet 1 Sg°^ and triplet states 3 S can be used 
as two qubit states in quantum computing. In a pro- 
posed quantum SWAP gate^ the gate operation time 
t ~ X/Jg-T, where Js-t being the singlet-triplet energy 
splitting. The singlet-triplet splitting of M Q/ j is shown in 
Fig. [7] on a semi-log plot. We see that it decay approxi- 
mately exponentially with the inter-dot distance. We also 
show in Fig. [7] the singlet-triplet splitting of the homo- 
QDM M 77 . We found that the Js-t of homo-QDM 
M 77 is slightly smaller than the Js-t of hetero-QDM 
M a 0, though the hopping energies of M a p and M 77 are 
almost identical. In the hetero-QDM case, the singlet 
wavefunction has more weight on the lower energy dot 
and therefore lowers the singlet energy and increases the 
singlet-triplet splitting. 



B. Double occupation of one of the dots in a QDM 

Double occupation means that two electrons occupy 
the same dot in a QDM. If the double occupation rate is 
high, the quantum gate operation may fail. The double 
occupation rate also reflects the localization properties of 
electrons in the QDM. If the double occupation rate is 
zero, each dot has one electron, whereas double occupa- 
tion rate of 1 means that two electrons are always local- 
ize on a single dot. When the double occupation rate is 
0.5, two electrons are delocalized between two dots. The 
double occupation can be conveniently analyzed in the 
dot-localized basis by transforming the CI equations to 
the dot-localized basis. 8 In the simplest case, wc consider 
only the "s" orbital for each dot, which give six config- 
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urations as follows, \e T ,e B ), \ej,,e B ), \e T ,e B ), \e T ,e B ), is^ 
l e Lj e s) an d l e T> e T)- The Hamiltonian in this basis set 



6t + es 



H = 



- Jte 
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2ct + Jtt / 



(9) 



Double occupation rate vs. d 
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FIG. 8: The double occupation rate of the ground-state sin- 
glet 1 Eg <1 ' vs. inter-dot distance for (a) hetero-QDM M a p and 
(b) homo-QDM M 77 . 



where t — txB- We ignored in Eq.© the off-diagonal 
Coulomb integrals, which are much smaller than the hop- 
ping t. 

The calculation of the matrix elements of Eq. ([9]) is 
described in Sec lHI Bl The two electrons can be either 
both on the top dots, or both on the bottom dots, or one 
on the top and the other on the bottom dots. We denote 
by {Ixfp' Xi> p>)} the configuration where one electron is 
on the Z-th orbital of the p dot with spin cr, and the 
other electron is on the i'-th orbital of the p' dot with 
spin cr'. Then the double occupation rate Qpp* in the 
many-particle state v is the probability of two electrons 



occupying the dot p— (T or B) at the same time, i.e., 

Q<$= £ p »(\xl P >x?>J)> ( 10 ) 



la.V '■ 



where P v iC) is the weight of the configuration C in the 
many-body wave functions of state v. The total prob- 
ability of two electrons being on the same dot is then 

Qtl = Q' 



TT 



Qftl for the f-th state. 



BB 

We plot Qtot, Qtt and Qbb of state for M a p in 



Fig. [8ja) and for M 77 in Fig. [8jb). We also performed 
calculations on a "symmetrized" model QDM M a / a r by 
setting e' T = e' B = (ey + es)/2 and J' TT = J' BB = 
(Jtt + Jbb)/2 of M a p in Eq. ©. M Q / Q / represents 
an ideal homo-QDM, without the asymmetry caused by 
strain, size and alloy composition effects. When compare 
the double occupation of the hetero- and homo-QDMs, 
we see that (i) For both types of QDMs, Q to t ~ 0.5 at 
d ~ 4.5 nm, meaning that two electrons are delocalized 
on two dots. For both QDMs, Qtot decays monotonically 
with the inter-dot distance, and at d ~ 10 nm, Qtot ~ 0, 
meaning that the two electrons are about each localized 
on one of the two dots. 

On the other hand, the double occupation of individual 
dot Qtt and Qbb differ substantially for homo-QDMs 
and hctcro-QDMs: 

(ii) For the homo-QDM M a t a i, Qbb = Qtt and de- 
cay monotonically with the inter-dot distances. Qbb 
and Qtt of M 77 have similar features, although Qbb 
is slightly different from Qtt due to the strain and alloy 
effects. This feature is also seen in the homo-QDM made 
of pure InAs/GaAs dots&2. In the hetero-QDMs M af3 , 
Qtt behaves very differently from Qbb because the ef- 
fective single-particle energy ex < es ■ Whereas Qbb de- 
cays monotonically with the inter-dot distance, Qtt has 
a maximum at d ~ 7 nm. The reason is that at d ~ 4.5 
nm, the hopping energy 2t is much larger than es — &t, 
therefore the electrons can overcome the energy barrier 
between the top and bottom dots and distribute evenly 
between two dots, leading to Qtt ~ Qbb- At larger d, 
2t <C es — ex, and the electrons would prefer to localize 
on the top dots, leading to Qtt 3> Qbb- Therefore, even 
when the total double occupation rate drops down, Qtt 
still increases and reaches the maximum at d=7 nm. For 
d >7 nm, Qtt decays as Qtot decays. 
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DOE of two-electron states 
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FIG. 9: (Color online) The degree of entanglement of two- 
electron states 1 'E g , 1 E„, 1 Eg and 3 S„, in (a) the hetero- 
QDM M a p, (b) the homo-QDM iW T7 and (c) the model "sym- 
metrized" homo-QDM M a i a i. 



(iii) The homo-QDMs M a > a i and M 77 have almost the 
same total double occupation, both smaller than that of 
the hetero-QDM M a p. The asymmetry between two dots 
increases the total double occupation. In an extreme 
case, where er <§C eg, the two electrons could always 
localize on the top dots, leading to Qtot=QTT=l- 



quantum correlation function for two electrons which has 
similar properties as the DOE. 24 However, the general- 
ization of this quantum correlation function to a system 
that has more than two single-particle levels is compli- 
cated. We proposed a DOE measure 9 - for indistinguish- 
able fermions using the Slater decomposition: 



,24,31 



as, 



S 



^ lo s2 



(11) 



where, Zi are Slater decomposition coefficients and 
zf=l. As shown in Ref. M, the DOE measure Eq. lTTTj) 
reduces to the usual Von Neumann entropy for distin- 
guishable particles when the two electrons are far from 
each other. In Rcfs. I25ll26l . a similar DOE measure was 
defined, which however due to a different normalization 
condition for Zi was used, does not reduce to the usual 
Von Neumann entropy even when the two electrons can 
be distinguished by their sites. 

The DOE of E states calculated from Eq. ^TTJ for the 
hetero-QDM M Q( g, the homo-QDM M 77 , and the model 
homo-QDM M a r a i are shown in Fig. [9ja), (b) and (c) 
respectively. All of the three QDMs have the following 
features: (i) <S( 1 Eg a ' > ) is close to zero (unentangled) at 
d r~j 4.5 nm, and close to unity (fully entangled) at d ~ 
10 nm. (ii) <S(qE) is almost unity (fully entangled) at 
all inter-dot distances. However, <S( 1 £g°' ) ) of the homo- 
QDM Af 77 (which is very close the S^E^) of M a > a >) is 

larger than S^'Sg ) of the hetero-QDM M a p, showing 
that the asymmetry in a QDM lowers the two-electron 
entanglement of the ground state singlet. 



In contrast to 5( 1 E^ aJ ) and S(g£), S^Hg') and 
5( 1 E„) are very sensitive to the asymmetry of the QDMs. 
In general, if the two dots have identical electronic struc- 
tures (e.g., in the simple Hubbard model), 5( 1 E^) = 
S( l E g a) ) and 5( 1 E U ) = 1, 9 as is illustrated in Fig. GJc) for 
M a ' a i. For M 77 , which is somehow asymmetric due to 

the strain and alloy effects, 5( 1 E^ ) is close to 5( 1 Eg') 
at small d, and drops down at large d, whereas for M a p, 



S{ x Y, g 0> ) is different from S{ l Y> g a) ) at all inter-dot dis- 
tances. The slight asymmetry in M 77 also causes iS( 1 E u ) 
to drop down at large d, similar to 5( £„) of M a p. 



V. ENTANGLEMENT 
A. Degree of entanglement for two electrons 

The degree of entanglement (DOE) is one of the most 
important quantities for successful quantum gate oper- 
ations. For distinguishable particles such as an elec- 
tron and a hole, the DOE can be calculated from 
the Von Neumann-entropy formulatio m 20 ' 21 ! 22 ' 23 How- 
ever, Von Neumann entropy formulation can not be 
used directly to calculate DOE for indistinguishable 
particles i 24 ' 25 i 26 ' 27 i 28 ' 29 i 30 Schliemann et al. proposed a 



B. Degree of entanglement vs double occupation 

Experimentally, it is very hard to measure the DOE 
of two electrons in the QDM directly, while it is rela- 
tively easy to measure the possibility of double occupa- 
tion. Therefore it would be useful to explore the rela- 
tion between DOE and the double occupation rate. The 
triplet states 3 E have negligible double occupation rate 
due to the Pauli exclusion principle. Here, we discuss 
the relation between DOE and double occupation rate 
for the ground state singlet 1 T, g a \ We consider the sim- 
plest case, where only "s" orbital in each dot is consid- 
ered. The ground state singlet 1 Ti g a ' > wavefunction can 
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Degree of entanglement vs. double occupation 




Double Occupation rate 



rate of the above ideal model in Fig. [10] in a black solid 
line. We also present in the same figure, the DOE of 
M a p, M 77 and M a r a i vs. double occupation rate. We 
found that the double occupation dependence of DOE 
for the homo-QDM M a i a i has perfect agreement with 
the analytical result, which is also true for M~ ll even 
though it has small asymmetry in the molecule due to 
the strain and alloy effects. We also checked the homo- 
QDM made of pure InAs/GaAs dots&2, and found the 

same double occupation dependence of DOE for the 1 Y> < f^ 
state, indicating this is a robust feature for homo-QDMs. 
However, the double occupation dependence of DOE for 
M a p deviates from the ideal case because dots a and j3 
arc different. 



FIG. 10: (Color online) Comparison of the degree of entan- 
glement vs. double occupation rate for hetero- and homo- 
QDMs. The black solid line represents the analytical results 
of homo-QDM, and the red dashed line represents the numer- 
ical results for homo-QDM M a i a i, and Af 77 , whereas the blue 
line represents the results for hetero-QDMs M a p. 



be generally written as, 

^EW) = c 1 \el r ,e B )+c 2 \e B , 4)+c 3 |e T , e l T )+c 4 \e B , e B ) , 

(12) 

and |ci| 2 + |c 2 | 2 + |c 3 | 2 + |c 4 | 2 = 1. Alternatively, we have 



where, 



*( 1 EW)=^ Wii |i)®|i) 



/ -c 3 -ci 

c 3 c 2 

-c 2 -C4 

\ci c 4 



(13) 



U) = 



and |») ,\j) = |4), K), \e B ), \e%) 
use Eq. [IT] to calculate the DOE, where 



(14) 



We can 



1/2(1 - V 1 - 4 ( c iC2 - c 3 c 4 ) 2 ) and z\ 



1/2(1 



yl — 4(cic 2 — c 3 c 4 ) 2 ) are the eigenvalues of lo^ uj. For 
a QDM with reflection symmetry, we have c\ = c 2 and 
c 3 = c 4 , and therefore z\ = 1/2(1 - - (1 - 4c 2 ) 2 ), 
and z\ = 1/2(1+ \j\ - (1 - 4c§) 2 ). Using the definition 
of double occupation rate, Q t ot = c§ + c 2 , we have 



z\ = 1/2(1 - y/1 - (1 - 2Q tot ) 2 ) , 

4 = 1/2(1 + v / l-(l-2Qtot) 2 ). (15) 

The DOE of 1 Sg Q ^ is calculated by substituting z 2 , z 2 
into Eq. (fTT|) . We plot the DOE vs. double occupation 



VI. SUMMARY 

We have studied the electronic structures of quantum 
dot molecules made of (In,Ga)As/GaAs dots of different 
sizes (hetero-QDMs), and compare them to that of quan- 
tum dot molecules made of identical dots (homo-QDMs). 
We found that while the hetero-QDMs and homo-QDMs 
have relatively similar electronic structures at short inter- 
dot distance, they differ significantly at large inter-dot 
distance, (i) Unlike those of homo-QDMs, the single- 
particle molecular orbitals of hetero-QDMs convert to 
dot localized orbitals at large inter-dot distance, (ii) Con- 
sequently, the bonding-antibonding splitting of molecu- 
lar orbitals is significantly larger than the electron hop- 
ping energy in a hetero-QMD at large inter-dot distance, 
whereas for homo-QDM, the bonding-antibonding split- 
ting is very similar to the hopping energy, (iii) The asym- 
metry of the QDM will significantly increase the dou- 
ble occupation for the two-electron ground states, and 
therefore lowers the degree of entanglement of the two 
electrons. 
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